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ABSTRACT 

The infrared (IR) emission of "M* galaxies" (10 10 4 < M stal < 10 M Q ) in galaxy pairs, derived using data 
obtained in Herschel (PEP/HerMES) and Spitzer (S-COSMOS) surveys, is compared to that of single disk 
galaxies in well matched control samples to study the cosmic evolution of the star-formation enhancement in- 
duced by galaxy-galaxy interaction. Both the mean IR SED and mean IR luminosity of star-forming galaxies 
(SFGs) in SFG+SFG (S+S) pairs in the redshift bin of 0.6 < z < 1 are consistent with no star-formation en- 
hancement. SFGs in S+S pairs in a lower redshift bin of 0.2 < z < 0.6 show marginal evidence for a weak 
star-formation enhancement. Together with the significant and strong sSFR enhancement shown by SFGs in a 
local sample of S+S pairs (obtained using previously published Spitzer observations), our results reveal a trend 
for the star-formation enhancement in S+S pairs to decrease with increasing redshift. Between z = and z = 1, 
this decline of interaction-induced star-formation enhancement occurs in parallel with the dramatic increase (by 
a factor of ~ 10) of the sSFR of single SFGs, both can be explained by the higher gas fraction in higher z disks. 
SFGs in mixed pairs (S+E pairs) do not show any significant star-formation enhancement at any redshift. The 
difference between SFGs in S+S pairs and in S+E pairs suggests a modulation of the sSFR by the inter-galactic 
medium (IGM) in the dark matter halos (DMH) hosting these pairs. 

Subject headings: galaxies: interactions — galaxies: evolution — galaxies: starburst — galaxies: general 



1. INTRODUCTION 

In some cosmological simulations (Guiderdoni et al. 
1998; Somerville et al. 2000; Baughetal. 2005), it is as- 
sumed that merger-induced star-formation is the major 
(or even the dominant) contribution to the high star- 
formation rate (SFR) at z ~ 1 - 2, the epoch when the 
SFR density in the universe peaks. Results of early 
HST surveys (Driver et al. 1995; Glazebrook et al. 1995; 
Abraham et al. 1996; Brinchmann et al. 1998; LeFevre et al. 
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2000; Conselice et al. 2003) were indeed consistent with this 
assumption. However, more recent observations (Bell et al. 
2005; Noeskeetal. 2007; Elbaz et al. 2007; Daddietal. 
2007; Lotz et al. 2008; Pannella et al. 2009; Rodighiero et al. 
2011; Pengetal. 2010; Elbaz et al. 2011; Wuyts et al. 
2011) have favored a scenario in which the so-called "main 
sequence" (MS) normal star-forming galaxies (SFG), contin- 
uously fueled by smooth accretion of cold intergalactic gas 
("cold flows", Dekel et al. 2009; Keres et al. 2009), dominate 
the cosmic SFR ever since z~ 4 while the contribution from 
merger-induced starbursts (the outliers in the SFR-M star plot) 
plays only a minor role (Rodighiero et al. 2011, but see 
Bridge et al. 2007; Shi et al. 2009). 

On the other hand, it is possible that much of the SFR 
in the MS galaxies could be associated with galaxy merg- 
ers. Many recent studies on merger rate evolution found 
that the fraction of galaxies in mergers increases significantly 
from z=0 to z=l (Kartaltepe et al. 2007; Conselice et al. 2009; 
Bridge et al. 2010; Xu et al. 2012, and ~ 10% galaxies with 
Z ~ 1 are in close major-merger pairs (see Lin et al. 2008; 
Lotz et al. 2008 for different results). In the local universe, 
the average SFR of merging galaxies (as found in opti- 
cal/NIR selected pairs) is only a factor of ^ 2 - 3 of that 
of single disk galaxies (Kennicutt et al. 1987; Ellison et al. 
2008; Xu et al. 2010). Therefore most of these galaxies be- 
long to the MS population, while extreme starbursts such as 
those in local ultra-luminous galaxies (ULIRGs) are very rare 
(Sanders & Mirabel 1996). If mergers of higher redshifts have 
similar or higher level of SFR enhancement, then ^ 30% of 

the SFR associated with MS galaxies could be due to mergers. 

In this work, exploiting the infrared (IR) data obtained 
using Herschel Space Observatory (Pilbratt et al. 2010) and 
Spitzer Space Telescope (Werner et al. 2004), we study the 
star-formation enhancement of SFGs in close major-merger 
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pairs since z=l. Throughout this paper, we adopt the 
A-cosmology with fl m = 0.3 and f^A = 0.7, and Ho = 
70 (km sec _1 Mpc _1 ). 

2. SAMPLES IN THE COSMOS FIELD 

Two samples of paired SFGs with 0.2 < z < 1 .0 are studied 
in this work, both are confined to "M* galaxies" with stellar 
mass in the range of 10 10 4 < M sta r < 10 11 M Q . According to 
Ilbert et al. (2010), the "M*" (the turning point in the galaxy 
stellar mass function) of SFGs is rather constant against red- 
shift since z = 2 : M* ~ 10 10 8 Mq for "intermediate active" 
galaxies and M* ~ 10 105 Mq for "high active" galaxies. The 
"S+S" sample includes 124 SFGs in 62 SFG+SFG pairs and 
the "S+E" sample 44 SFGs in 44 mixed pairs, all having red- 
shifts in the range of [0.2,1.0]. 

Galaxies in both samples are taken from the sample of close 

(5 < r pro j < 20 h" 1 kpc)major-mergerpairs (M^/M^ < 2.5) 
in the COSMOS field (CPAIR, Xu et al. 2012). The par- 
ent sample of CPAIR was selected from the photo-z catalog 
of the COSMOS field (Capak et al. 2007; Ilbert et al. 2009; 
Drory et al. 2009) with stellar mass limits of log(M min /M Q ) = 
[9.0,9.4,9.8,10.2] in four redshift bins that equally divide 
the redshift range of [0.2,1.0], respectively. These limits 
are above the completeness limits for the stellar mass of 
the photo-z sample (Drory et al. 2009). It was required in 
the CPAIR selection that the photo-z's of the two galaxies 
in a pair candidate satisfy the following criterion: |Zp" ot — 

< 0-03, where z% )t and z% t are the photo- 
z's of the primary and secondary, respectively. The redshifts 
of both galaxies in a pair are then assumed to be the same, es- 
timated by the mean of the two photo-z's. The classification 
of SFGs (S) or "passive galaxies" (E) was based on the SED 
fitting and taken from Drory et al. (2009). 

Both samples are divided into 2 redshift bins: [0.2,0.6] 
(low-z bin) and [0.6,1.0] (high-z bin). These bins are 2 times 
wider than those in the original CPAIR sample. Making 
the size of subsamples in redshift bins larger allows better 
statistics. Since for most pairs the beams of IR observa- 
tions (FWHM > 5") cannot resolve a pair into two component 
galaxies, the S+S sample selection was carried out on pairs in- 
stead of on individual galaxies: a pair will be included in the 

sample if [log(M^/M )+log(M^/M )]/2 is in the range 
of [10.4, 11]. The S+E sample includes only the SFGs in the 
mixed pairs and excludes the early-type galaxies. There are 7 
(11) S+S pairs (S+E pairs) in the low-z bin with the median 
redshift of 0.45 (0.41), and 55 (33) in the high-z bin with the 
median redshift of 0.81 (0.82). 

The two control samples, one for each sample of paired 
SFGs, are selected from the same parent sample of CPAIR. 
In order to minimize the statistical errors in our final results 
that are due to the controls, each paired galaxy was matched 
by 10 control galaxies (none being included more than once). 
The following criteria were used in the selection of the con- 
trols: (1) redshift match: |z^ t - zj*J/(l + z|| ot ) < 0.03, 

where z^,"' is the photo-z of the control candidate, and z P p hot 
the photo-z of the paired galaxy to be matched; (2) mass 
match: | logtM!;™ 1 ) - log(MP g al .) | < 0.1, where M^ l t is the stel- 
lar mass of the control candidate, and M^ ot the stellar mass 
of the paired galaxy to be matched; (3) type match: all con- 
trol galaxies are SFGs; (4) local density match: |p££ nt — pJJbl < 

MAX(2, Jtf^), where 

Pnb" i s th e neighbor counts (for neigh- 



bor galaxies of log(M star /MQ) > 10.2) within IMpc (comov- 
ing) of the control candidate, and p^ the neighbor counts 
around the paired galaxy. Here, in order to have adequate 
numbers of control candidates for paired galaxies with very 
low local densities, we set a minimum (min = 2) to the den- 
sity match limit which is otherwise equal to the Poisson error 

In order to exclude interacting galaxies, the following two 
additional criteria were applied to the selection of control 
candidates: (5) no companion within the projected distance 
r p = SOh-'kpc that has |<5z phot |/(l + z£°£) < 0.03 and more 
massive than 0.2 x M™^; (6) A < 0.35 and G+0.4 x A < 0.66, 
where A is the asymmetry parameter (Conselice et al. 2000) 
and G the Gini Coefficient (Lotz et al. 2008). The criterion (6) 
excludes morphologically disturbed galaxies that are gener- 
ally associated with mergers (Lotz et al. 2008; Conselice et al. 
2009). Visual inspections of 100 galaxies randomly picked 
from the control samples showed no obvious merger candi- 
dates except for a few ambiguous cases of galaxies with pos- 
sible weak and faint distortion features. 

3. LOCAL SAMPLES 

In order to provide the local benchmarks for the evolution- 
ary study, two local samples of paired SFGs, one for those in 
S+S pairs and another for those in S+E pairs, and their corre- 
sponding control samples are also studied. These samples are 
adopted from Xuet al. (2010), who observed a nearly com- 
plete K-band selected pair sample (KPAIR) using Spitzer in 
the 3.6, 4.5, 5.8, 8, 24, 70, and 160 pm bands. The sample 
was selected from cross-matches between 2MASS and SDSS- 
DR3 galaxies, including 27 S+S and S+E pairs that have 

5 < r proj < 20 h" 1 kpc and mass ratio of M^/M^ < 2.5. The 
selection criteria are nearly the same to those for the CPAIR 
sample. The paired galaxies are then one-to-one matched (ac- 
cording to the stellar mass and redshift) by a sample of con- 
trol galaxies selected from the SWIRE survey (Lonsdale et al. 
2003) and the SINGS survey (Kennicutt et al. 2003). There 
are 39 non-AGN SFGs in the pair sample of Xu et al. (2010). 
Among them, all paired SFGs (and their controls) in the mass 
range of 10 10 ' 4 < M star < 10 110 M Q (stellar mass estimated 
from the K-band luminosity using Kroupa IMF) are included 
in the local samples of this study. There are 20 (7) SFGs in 
the local S+S (S+E) sample, and the same number of galaxies 
in the corresponding control sample. 

4. THEIR DATA 

The COSMOS field was surveyed by the PACS Evolution- 
ary Probe (PEP, Lutz et al. 201 1) using PACS photometer ar- 
rays (Poglitsch et al. 2010). The 3-cr sensitivities in the PACS 
100 and 160 pm maps are 5.0 mJy and 10.2 mJy, respec- 
tively (Berta et al. 201 1). The 250, 350 and 500 pm observa- 
tions were carried out by Herschel Multi-tiered Extragalac- 
tic Survey^ (HerMES, Oliver et al. 2012) using SPIRE pho- 
tometer arrays (Griffin et al. 2010; Swinyard et al. 2010). The 
maps in the three SPIRE bands (Levenson et al. 2010) are 
confusion limited at 3-cr level of 17.4, 18.8, and 20.4 mJy, 
respectively (Nguyen et al. 2010). In addition, we also in- 
cluded the Spitzer 24 pm data from the S-COSMOS survey 
(Sanders et al. 2007). The 24 pm source catalog of the COS- 
MOS field (LeFloc'h et al. 2009) is flux limited at f Mfim = 
80 piy. 

14 http://hermes.sussex.ac.uk 
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Table 1 

Characteristics of the IR Data 
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4/11 
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Note. — a) Insrument; b) full width at half maximum of the point spread function (PSF); c) flux limit; d) signal-to-noise ratio at the flux limit; e) low-z bin 
(0.2 < z < 0.6); /) high-z bin (0.6 < z < 1). 



We assume that the IR emission in an S+E pair is due to the 
SFG component only, ignoring the contribution from the E 
component (Domingue et al. 2003). Because most S+S pairs 
(median angular separation ~ 3") are not resolved in the IR 
bands, each S+S pair is treated as a single source. Table 1 
provides the main characteristics of the IR data. In general, 
the detection rates of our samples in the 5 Herschel bands are 
rather low. In particular, except for S+S pairs, sources in the 
high-z bin (0.6 < z < 1) have detection rates < 10% in all 
Herschel samples. The S+S pairs have slightly higher detec- 
tion rates because most of them are not resolved by Herschel, 
therefore each Herschel source includes contributions from 
both component galaxies. The detection rates in the MIPS 24 
fj,m band are between 45.5% and 85.7%. 

5. STACKING ANALYSES AND RESULTS 

Because of the low IR detection rates, we studied the IR 
emission of our sources through stacking. Two stacking meth- 
ods are considered. The first is a modified version of "cleaned 
stacking". In its original form (Zheng et al. 2007), there are 
four steps in this method: (1) detect and extract all bright 
sources above a given detection threshold from the observa- 
tional image; (2) separate the target sources into detected and 
undetected by matching them with the catalog of detected 
sources; (3) stack the images of undetected sources, which 
are cut from the residual map with the detected sources ex- 
tracted, and measure the flux; (4) derive the mean flux of 
target sources f = (X)i f det,i + f stack) /Ntotai, where f det .i is the 
flux of the i-th detected source, f stac k the flux measured in the 
stacked image, and Ntotai the total number of sources in the 
target sample. Compared to uncleaned stacking (Dole et al. 
2006), cleaned stacking suffers less confusion problems due 
to unrelated background sources. However, the result might 
be biased due to flux boosting for detected sources, in partic- 
ular if the detection limit is set at a low s/a level. In order to 
avoid this bias, we made the following modifications. First, 
the detection thresholds were set relatively high (at s/a ~ 6) 
when doing the detection and extraction on observed images 
using StarFinder (Diolaiti et al. 2000). Then, the extractions 
corresponding to detected target sources were restored back 
to the residual map; the subsequent stacking, using the resid- 
ual map with detected background sources extracted, went 
through all target sources including both detected and unde- 
tected. Finally, the mean flux of a given band for a given 
sample was derived by f = f s tack/Nt tai- 

Some sources in our samples are marginally resolved in the 
IR images, in particular the S+S pairs with sep > 5" (~ 7% 
of all S+S pairs). In order to ensure that a stacked image in- 
cludes all fluxes from all sources, we searched in the catalog 



of detected sources for all matches to a given S+S pair using 
a searching radius of 0.5 x (FWHM 2 + sep 2 ). For SFGs in 
the S+E sample and two control samples, which are generally 
point sources in all IR bands (e.g. the mean Petrosian radius 
of galaxies in the S+S control sample is 1 .2" with a disper- 
sion of 0.6"), the searching radius is 0.5 x FWHM. In the 
photometry of the stacked images, we choose to use a con- 
stant aperture of radius R = 1 1 .2" for the three MIPS and PEP 
bands at 24, 100 and 160 \im. The corresponding aperture 
corrections are 1.09, 1.05, 1.21, respectively, according to our 
empirical tests using StarFinder. 

In the second method, the "covariance" method 
(Marsden et al. 2009), the mean flux of a sample of tar- 
get sources is estimated by f = X^k^k/N, wnere Dk is the 
measured flux density (in units of Jy/beam) in the map pixel 
that contains the k-th source in the target sample. The map, 
without any cleaning of unrelated sources, is background 
subtracted with zero mean. The method explicitly assumes 
that target sources are randomly distributed in the sky 
according to the Poisson statistics, which is certainly not 
valid for paired galaxies, therefore it should not be applied 
to individual galaxies in S+S pairs. Also, this method may 
significantly under-estimate the mean flux of target sources if 
some of them are resolved (even only partially). 

Our experiments showed that for SPIRE images at 250, 
350 and 500 /jjn, which are predominantly confusion limited 
(Oliver et al. 2012), the mean fluxes derived using the covari- 
ance method and cleaned stacking are all consistent with each 
other within 1<t error; and results of the covariance method 
generally have slightly lower errors than those of cleaned 
stacking. On the other hand, in the MIPS 24 pm and PACS 
100 fim bands, the mean fluxes of S+S pairs derived using 
the covariance method are systematically lower than those ob- 
tained using cleaned stacking, presumably due to missing flux 
of resolved sources in the former. Therefore, in the final anal- 
yses, we chose to use the cleaned stacking method for the 24, 
100 and 160 fjm bands, and the covariance method for the 
250, 350 and 500 fim bands. 

The results are listed in Table 2. All errors were estimated 
by bootstrapping (Efron 1979). For a given sample and band, 
the error includes both the statistical dispersion of the sam- 
ple and the measurement noise (instrumental and confusion 
noise). A 3-cr upper-limit is assigned to a mean flux if it is 
less than fr lm / \/N, where fr, m is the 3-cr sensitivity limit of the 
band, and N the size of the sample. For the S+S sample, N 
is the number of pairs; for other samples, N is the number of 
galaxies. For SFGs in the S+S sample, the mean fluxes, errors, 
and upper-limits were all obtained by dividing corresponding 
values of the pairs by a factor of 2. 
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Table 2 

IR Emission of Paired Galaxies and Control Galaxies a 
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log(LiR) 
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7 


0.45 
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9.78 ±3.99 
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< 3.86 c 
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55 


0.81 


0.13 ±0.02 


2.43 ±0.66 


5.00 ±0.94 


5.50±0.82 


4.56 ±0.79 


3.23 ±0.60 


11.15±0.06 


0.04 ±0.07 


Control 


140 


0.45 


0.28 ±0.03 


5.75 ±0.61 


10.25 ±0.91 


8.75 ±0.73 


5.71 ±0.71 


2.58 ±0.64 


10.86 ±0.04 
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0.82 


0.13 ±0.01 


2.03±0.19 


4.13 ±0.24 


5.16±0.27 


4.32 ±0.23 


2.34±0.21 


11.11 ±0.03 




S+E 


11 


0.41 


0.35 ±0.08 


8.69 ±2.63 


13.01 ±4.47 


13.18± 3.74 


8.43 ±2.93 


<6.15 c 


10.96±0.11 


0.05 ±0.11 




33 


0.82 


0.14 ±0.03 


1.76 ±1.07 


4.35 ±1.85 


5.95 ±1.65 


< 3.27 c 


< 3.55 


11.13±0.12 


0.04±0.12 


Control 


110 


0.41 


0.39 ±0.04 


8. 19 ± 1.07 


14.06± 1.38 


11.45 ± 1.02 


6.18±0.75 


2.51 ±0.69 


10.90 ±0.04 






330 


0.82 


0.14±0.01 


1.95 ±0.22 


3.40 ±0.37 


5.13 ±0.43 


4.26 ±0.46 


2.66 ±0.41 


11. 09 ±0.04 





Note. — a) Mean IR fluxes were estimated through stacking, and the errors through bootstrapping; b) mean sSFR enhancement index; c) 3-cr upper-limit (= fii m /\/N, 
where f|j m is the 3-<r sensitivity limit of the band, and N the size of the sample); d) mean IR fluxes, errors and upper-limits for paired SFGs were derived by dividing 
corresponding values of S+S pairs by 2. 
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Figure 1. Left Panel: The mean IR SED of SFGs in the S+S sample and of those in the S+S control sample. The red curve is the best fit to data of SFGs in 
S+S pairs with 0.2 < z < 0.6, using the semi-empirical SED of Arp 244 (the Antennae Galaxies) taken from Xu et al. (2001). The other three curves are the 
best fits to data of SFGs in S+S pairs with 0.6 < z. < 1.0 and SFGs in the S+S control sample with 0.2 <z. < 0.6 and 0.6 < z. < 1.0, respectively, using the same 
semi-empirical SED of normal Sc galaxy NGC 6181 (Xu et al. 2001). Right Panel: The mean IR SED of SFGs in the S+E sample and of those in the S+E 
control sample. All curves are the best fits to data using the same semi-empirical SED of normal Sc galaxy NGC 6181. 



We neglected the bias due to clustering of IR sources 
(Bethermin et al. 2012). Our experiments using simulations 
based on the algorithm described in Bethermin et al. (2012) 
showed that for our (mass-selected) samples the clustering 
bias is not strong, ~ 15% for the mean fluxes in the 500 fim 
band and ~ 7% for those in other bands, without significant 
dependence on redshift. Because of the local density match 
between the control samples and the pair samples (Section 
2), the biases for pairs and for controls should be the same 
and therefore cancel each other out in the star-formation en- 
hancement analysis. On the same ground, any other possible 
systematic biases in the mean fluxes (e.g. due to filtering of 
the maps) have also been neglected. 

6. STAR-FORMATION ENHANCEMENT IN PAIRED SFGS 

Using the median redshift of each subsample, the mean IR 
fluxes obtained by stacking are converted to the rest-frame 
vh y and plotted in Fig. 1. The semi-empirical SEDs in the 
SED library of Xu et al. (2001) were searched for the best fits 
to these data. The SED library 15 is based on the IRAS, ISO, 
and ground based sub-mm observations of 837 IRAS 25 /ira 
band selected galaxies. Among bright FIR sources (fioo/im > 
10 Jy) in the library, we found that the SED of local merger 
Arp 244 (the Antennae Galaxies) fits best the data of SFGs in 

15 http://spider.ipac.caltech.edu/staff/cxu/sed_lib/intro_sed_lib.html 



S+S sample in the low-z bin (0.2 < z < 0.6). Interestingly, the 
mean SEDs of all other subsamples, including the high-z bin 
(0.6 < z < 1) of the S+S sample and both redshift bins of the 
S+E sample, have similar shapes and can be best fitted by the 
SED of the same local normal Sc galaxy NGC 6181 (Fig. 1). 

The inverse variance weighted least-squares fit to data in 
the 100, 160, 250, 350 pm bands (i.e. data close to the IR 
peak), using the corresponding best-fit SED, was carried out 
for each subsample to determine the mean total IR luminosity 
Lir(5 - 1000/im). The variance of Lir was estimated by 

where all summations are over the four bands at 100, 160, 
250, 350 fim, and Covj.j is the covariance between band i and 
band j. If the errors in different bands were independent with 
each other, i.e. Covj j = for i j, the variance of Lir could 
have been estimated by a 2 = 1/of. However for our 

data we found rather significant covariance between different 
bands (also estimated using bootstrapping). 

Resulting Lir's are listed in Table 2 and plotted in Fig. 2 
against redshift. Also plotted are mean IR luminosities of 
SFGs in the local samples. The mean Lir's of the local sam- 
ples are derived using Spitzer data at 8, 24, 70, and 160 fim 
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Figure 2. Plot of mean log(LiR) vs. redshift. Each data point and the error 
bar at z > were derived through fitting the corresponding SED in Fig. 1. 
Data points for the local (z = 0) samples were obtained using Spitzer observa- 
tions of Xu et al. (2010). There are 20 and 7 SFGs in the local S+S and S+E 
samples. Their control samples have the same sizes. 
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Figure 3. Plot of mean sSFR enhancement index, e = log(sSFR pll j, ) — 
log(sSFRcom). vs. redshift. Because each pair sample and its control sam- 
ple have the same mean stellar mass, e = log(L[R pll i r )-log(LiR cont ). 

(Xu et al. 2010). Fig. 3 shows the sSFR enhancement of SFGs 
in S+S pairs and in S+E pairs. The sSFR enhancement index 
is defined as 



e = log(sSFR pair ) - log(sSFR a 
= log(L IRiPair ) - \og{L IRtCont ). 



,) 



controls are both best-fitted by the same SED of a normal disk 
galaxy (NGC 6181). In the low-z bin (0.2 < z < 0.6), mainly 
because of the small number of S+S pairs in the bin, the re- 
sults are less clear-cut: their average sSFR enhancement index 
is e = 0.20 ±0.13, indicating a weak enhancement at rather 
low significance level of 87% (i.e. at 1.5-er level). The mean 
SED is best fitted by that of a galaxy merger (Arp 244) but, 
given the large error bars of individual data points, a fit with 
the SED of the normal spiral NGC 6181 cannot be ruled out. 
For SFGs in S+E pairs in both the low-z and high-z bins, re- 
sults on both the sSFR enhancement and the SED comparison 
are consistent with no enhancement in their star-formation ac- 
tivities. 

7. DISCUSSION 

7.1. S+S Pairs: Less Star-Formation Enhancement at Higher 

z 

The non-enhancement of the star-formation in SFGs in S+S 
pairs with z ~ 0.8 and the rather weak star-formation enhance- 
ment in those with z ~ 0.4, in contrast with the strong and 
significant star-formation enhancement (e = 0.50 ±0.16) in 
the local sample (Fig. 3), reveal a clear trend for the star- 
formation enhancement to decrease with redshift. It is un- 
likely that this trend is due to an increase of the spurious pairs 
fraction (SPF) with redshift. The probability for selecting (un- 
physical) projected pairs is significantly reduced by CPAIR 

selection criteria 5 < r proj < 20 h" 1 kpc and M^/M 2 ^ < 2.5, 
confining the companion search to a very small sky area 
(~ 10 arcsce 2 ) and a narrow mass range. Indeed, as shown 
by the Monte Carlo simulations (Xu et al. 2012), the spuri- 
ous pair fraction due to projected pairs (given the photo-z er- 
ror of the COSMOS survey and the pair selection criterion 
Sz p hot/(l +z p hot) < 0.03) is only 7 ± 3% at z=0.3 and 9 ± 3% 
at z = 0.9. Even after including the effect of galaxy clustering 
and counting pairs with <5V > 500 km s" 1 as spurious pairs, 
the spurious pair fraction is still rather low in the CPAIR sam- 
ple: SPF = 0.20 ± 0.06 at z = 0.3 and SPF = 0.22 ± 0.06 at 
z = 0.9. The local (z = 0) pairs were selected based on spectro- 
scopic redshifts (Xu et al. 2010). For them, SPF = 0.06±0.05 
(Xu et al. 2012), due to unbound galaxies that are clustered to- 
gether (therefore having recession velocity difference as low 
as 5V < 500 km s" 1 ). Statistically, the effect of the contami- 
nation of spurious pairs can be corrected as following 

(4) 



(2) 
(3) 



The last equation is valid because, by design, each pair sample 
and its control sample have the same mean stellar mass. Also, 
we neglected the uncertainties of the Lir as an SFR indicator 
that are due either to the dust heating by old stars or to the UV 
radiation of young stars escaping the dust absorption. These 
uncertainties should not affect the calculation of e because the 
effects on the paired SFGs and on the controls should be the 
same and therefore should cancel each other out. 

For SFGs in S+S pairs, data are consistent with un- 
enhanced star-formation in paired galaxies in the high-z bin 
(0.6 < z < 1). They show an average sSFR enhancement in- 
dex of e = 0.04 ± 0.07, and their mean SED and that of their 



e true = e obs /(l-SPF). 



For SFGs in S+S pairs in the high-z bin (0.6 < z < 1), this 
correction results in a small change in the mean enhancement 
index, from e Q b s = 0.04 to e true = 0.05. In summary, the low 
spurious pair fractions in our pair samples shall not affect the 
results on the sSFR enhancement significantly. 

In parallel to the decline of the sSFR enhancement in S+S 
pairs, mean Lir of normal SFGs in the control sample in- 
creases rapidly with redshift (Fig. 2). Strong SFR evolu- 
tion in normal disk galaxies since z ~ 1 has been well docu- 
mented (Bell et al. 2005; Zheng et al. 2007; Elbaz et al. 2007; 
Noeske et al. 2007), and has been attributed to the increase 
of gas fraction with redshift in these galaxies. Indeed, CO 
observations of massive SFGs at z ~ 1 show significantly 
higher molecular gas content compared to local disk galax- 
ies (Tacconi et al. 2010). According to Hopkins et al. (2009), 
the ratio between mass of gas consumed by an interaction- 
induced nuclear starburst and the disk gas mass decreases with 
increasing disk gas fraction as fburst/fgas oc (l-f gas ), where 
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fburst is the gas fraction (among the total baryonic mass) con- 
sumed by the starburst and f gas the disk gas fraction. This is 
because the gravitational torque imposed by the stellar disk to 
the gas disk is less effective when f gas is high, therefore less 
disk gas can sink to the nuclear region by losing angular mo- 
mentum to stars. We argue that the negative cosmic evolution 
of the sSFR enhancement in S+S pairs is due to the increase 
of f gas with redshift. 

The sSFR enhancement in local (z = 0) close major- 
merger pairs has been well established (Kennicutt et al. 1987; 
Xu & Sulentic 1991; Barton et al. 2000; Nikolic et al. 2004; 
Ellison et al. 2008). Meanwhile, Xu et al. (2010) found that 
the enhancement only occurs in massive SFGs ( ^ 10 10 M Q ) 

in S+S pairs. Paired SFGs with low M star have in gen- 
eral higher f gas than massive SFGs and show no signifi- 
cant sSFR enhancement. It is worth noting that studies of 
B-band selected interacting galaxies, which are biased for 
low mass SFGs, found no or very weak sSFR enhancement 
(Bergvall et al. 2003; Knapen & James 2009). It is plausible 
that the high f aas ( £ 30%) plays a significant role in the non- 
enhancement of the sSFR in low mass paired SFGs. 

There have only been a few previous statistical studies of 
the SFR of interacting galaxies in the universe at z > 0.1. 
Using UV and optical SFR indicators, which are prone to 
dust attenuation, both LeFevre et al. (2000) and Bridge et al. 
(2010) found that the SFR of interacting galaxies at z rS 1 are 

enhanced by a factor of ~ 2 compared to the field galax- 
ies, and the SFR enhancement increases with the redshift. 
On the other hand, LeFevre et al. (2000) noticed a difference 
between "upcoming major-mergers" and "ongoing major- 
mergers"; the former do not show significant SFR enhance- 
ment (as measured by [Oil] equivalent width). Bridge et al. 
(2010) also found that, in their sample of morphologically 
selected mergers, early stage mergers (close pairs with tidal 
bridges) have SFRs similar to field galaxies. These are consis- 
tent with our results of the non-enhancement for paired SFG at 
Z = 0.8, because our pair selection favors early stage mergers. 
However, different from early stage mergers in LeFevre et al. 
(2000) and Bridge et al. (2010), SFGs in S+S pairs at z = 
do show significant sSFR enhancement in this study. Given 
that our pairs at different redshifts are selected using the same 
selection criteria, it is unlikely that the ratio between the 
early stage and late stage mergers changes with the redshift. 
Lin et al. (2007), comparing the median LiR/M stal - (Lir being 
derived from the 24/im flux) of a sample of close pairs and 
that of pseudo-pairs constructed using conrol galaxies, found 
an enhancement of a factor of ~ 2 for pairs in the redshift 
range 0.1 < z < 1 . However, by assigning a very low Lir value 
(1O 9 L ) to undetected control galaxies, Lin et al. (2007) may 
have under-estimated the median L IR /M star of control pairs 
and over-estimated the enhancement of the close pairs (most 
being unresolved in the 24/im survey). Our result on the mean 
24/1/7Z emission of S+S pairs in the high-z bin differs from 
that of Lin et al. (2007), showing no enhancement compared 
to the controls (Table 2). Jogee et al. (2009) studied the SFR 
of morphologically selected mergers with 0.24 < z < 0.80, 
using both the UV and MIR (24 fjm) data, and found modest 
enhancement compared to normal galaxies. Studying an IR 
selected sample observed by Herschel, Hwang et al. (2011) 
found a factor of 1.8 - 4.0 sSFR enhancement for galaxies in 
close SFG-SFG pairs with redshifts between and 1 .2, with 
little variation in the enhancement versus redshift. However, 



given the blending effect, the true enhancement indicated by 
results of Hwang et al. (201 1) is significantly less (by a factor 
of ~ 0.5). Also their results are consistent with a decreas- 
ing star formation enhancement at higher z because higher z 
sources are likely to be affected more severely by blending. 

7.2. Cosmic Evolution of Extreme Starbursts 

According to simulations (Scudderet al. 2012), samples of 
close pairs (5 < r pro j < 20 h" 1 kpc) include mostly merging 
galaxies undergoing the first close encounter and those just 
before the final coalescence, both stages lasting ~ 10 8 years. 
They miss final stage mergers (with the time scale of ~ 10 7 
years) that are too close to be identified as binaries, and those 
already coalesced. Early IRAS observations (Kennicutt et al. 
1987; Sanders & Mirabel 1996) already revealed that final 
stage mergers show much stronger star-formation enhance- 
ment than paired galaxies, and extreme starbursts such as 
ULIRGs occur almost exclusively within these final stage 
mergers. Given the selection effects for pair samples, the cos- 
mic evolution of these extreme starbursts is not probed by the 
data in this study. 

On the other hand, there are indications in the literature that 
properties of extreme starbursts in final stage mergers might 
also change with redshift. Rujopakarn et al. (2011) found 
that ULIRGs at z ~ 1 and sub-millimeter galaxies (SMGs) at 
z ~ 2 have much more extended star-formation distributions 
than local ULIRGs, and argued that these are isolated galax- 
ies with high Esfr- However, studying the morphology of 70 
fj,m sources in the S-COSMOS survey, Kartaltepe et al. (2010) 
concluded that major mergers dominate the ULIRG popula- 
tion at z 1 . The CO and radio continuum observations 

of Biggs & Ivison (2008) and Engel et al. (2010) demonstrate 
that SMGs are also dominated by major mergers even though, 
different from local ULIRGs that typically have sizes of ;$ 1 

kpc (Rujopakarn et al. 2011), the extreme starbursts in high- 
z mergers are much more extended (~ 5 kpc). These au- 
thors suggested that the difference between the size of local 
ULIRGs and that of high-z SMGs is due to the selection ef- 
fect in the sense that SMGs are preferentially earlier stage 
mergers than local ULIRGs. However, it is possible that the 
star-formation in high-z SMGs (and high-z ULIRGs) is less 
centrally peaked because the strength of the nuclear starburst 
is suppressed due to the higher gas fraction in the disk, accord- 
ing to the same physical mechanism (Hopkins et al. 2009) 
that we invoked in the interpretation for the weakening of the 
sSFR enhancement in higher z S+S pairs. It will be interesting 
to check in future studies whether the percentage of extreme 
starbursts (with sSFR enhancement £ a factor of 100) among 

major mergers also decreases with redshift. 

7.3. SFGs in S+E Pairs: No sSFR Enhancement 

For SFGs in S+E pairs, the sSFR enhancement index e is 
consistent with at all redshifts. This agrees with the Spitzer 
results of Xuetal. (2010) which showed no sSFR enhance- 
ment for SFGs in S+E pairs in any stellar mass bin in a 
local close major-merger pair sample. Hwang et al. (2011) 
also found that SFGs with close early-type neighbors are not 
sSFR enhanced. The difference between SFGs in S+S pairs 
and in S+E pairs indicates that, in addition to gravitational 
tidal effects, the sSFR in a paired galaxy is influenced by 
the immediate surrounding environment. This hypothesis is 
in agreement with the correlation between sSFR's of the pri- 
maries and secondaries in major-merger S+S pairs (i.e. the 
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"Holmberg effect", Kennicutt et al. 1987; Xu et al. 2010). On 
the other hand, Xu et al. (2010) found no significant differ- 
ence between the local densities around S+S and S+E pairs 
in their local pair sample, which is confirmed again here for 
pairs with z = 0.2 - 1.0; the average counts of neighbors of 
M star > 1O 1O 2 M0 within 1 Mpc projected distance and with 
bzphot /(l +z p hot) < 0.007 (the la error of photo-z, Ilbert et al. 
2009) are 3.69 ± 0.45 and 3.73 ± 0.52 for S+S pairs and S+E 
pairs, respectively. Therefore, the linear scale of the environ- 
ment effect must be less than 1 Mpc. Because most galaxies 
in close major-merger pairs have entered the virial radius of 
the companions a long time ago, the two galaxies in a pair 
are likely to share the same IGM gas in a common dark mat- 
ter halo (DMH). The non-enhancement of the sSFR in SFGs 
in S+E pairs and the "Holmberg effect" for S+S pairs sug- 
gest a significant role of the IGM within a DMH to the sSFR 
of galaxies residing in the DMH. For example, when a DMH 
has strong (weak) "cold flows" (Dekel et al. 2009; Keres et al. 
2009), galaxies inside it may have abundant (scarce) cold gas 
supply to fuel active star-formation. This hypotheses seems 
particularly attractive because it can explain both the "Holm- 
berg effect" for S+S pairs and the non-enhancement of the 
SFG in S+E pairs. Another scenario for the IGM modula- 
tion of the sSFR involves the hot IGM gas which can strip 
the cold ISM gas in embedded galaxies (Park & Choi 2005), 
though this may have difficulty in explaining the "Holm- 
berg effect" for S+S pairs because detections of the hot IGM 
gas around spiral galaxies are very rare (Bensen et al. 2000; 
Anderson & Bregman 201 1). 
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